Spin polarons and molecules in strongly-interacting atomic Fermi gases 
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We examine pairing and molecule formation in strongly-interacting Fermi gases, and we discuss 
how radio-frequency (RF) spectroscopy can reveal these features. For the balanced case, the emer- 
gence of stable molecules in the BEC regime results in a two-peak structure in the RF spectrum 
with clearly visible medium effects on the low-energy part of the molecular wavefunction. For the 
highly-imbalanced case, we show the existence of a well-defined quasiparticle (a spin polaron) on 
both sides of the Feshbach resonance, we evaluate its lifetime, and we illustrate how its energy may 
be measured by RF spectroscopy. 
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The experimental realization of strongly-interacting 
Fermi gases with Feshbach resonances constitutes one 
of the major recent breakthroughs in ultracold atomic 
physics. As a result, we can now study the properties of 
Fermi gases as a function of interaction strength, tem- 
perature, and population imbalance (polarization). A 
fundamental question is how we should understand pair- 
ing in such strongly-interacting systems. In the limit 
where attractions are weak and attractive, the BCS in- 
stability can lead to the formation of Cooper pairs, which 
are at the origin of our understanding of superconductiv- 
ity. In the opposite Bose-Einstein condensation (BEC) 
limit, pairing leads to tightly bound molecules mixed 
with free atoms, which may be present due to either a 
nonzero temperature or polarization. In the intermedi- 
ate regime close to the unitarity limit, the many-body 
physics is much more complicated. Radio- frequency (RF) 
spectroscopy represents at present the main experimen- 
tal probe of strong correlations in ultracold atomic Fermi 
gases [l|, Q • Recent experiments performed on a gas in 
the superfluid phase seem to be well described by a two- 
body theory with the molecule energy as a fitting pa- 
rameter Q. However, microscopic calculations for the 
unpolarized state indicate that both the superfluid and 
normal states at unitarity cannot be described simply as 
molecules mixed with free atoms [H, Furthermore, the 
highly polarized state remains normal even at T = and 
presumably is a Fermi liquid Q. 

In this paper we investigate the nature of the strong 
correlations, pairing, and molecule formation for polar- 
ized and unpolarized atomic Fermi gases on both sides of 
the Feshbach resonance. Our analysis concerns the nor- 
mal phase, which arguably is not as well understood as 
the T — superfluid phase 0], and we analyze how sev- 
eral nontrivial effects can be revealed in RF spectroscopy. 
Pair fluctuations, which in the BCS and unitarity regimes 
lead to a pseudogap, change character in the BEC regime 
with the emergence of stable molecular states. We show 
how this change may be directly detected in the RF spec- 
trum. For the highly-polarized case, we demonstrate the 



existence of a well-defined quasiparticle (spin-polaron) in 
the minority spectral function on both sides of the Fesh- 
bach resonance. We calculate its lifetime, and show how 
its energy may be measured via RF spectroscopy. We 
further discuss how many-body effects significantly affect 
the molecular wavefunction and how this can be detected 
in RF spectroscopy. For a polarized gas, Pauli-blocking 
effects are shown to lead to an asymmetry between the 
minority and majority RF spectra. By comparing our re- 
sults to a simple model of atoms mixed with molecules, 
we present clear physical interpretations of our results. 
Our analysis is based on a microscopic theory with no 
fitting parameters and, as we show, for high imbalance it 
compares excellently with recent Monte-Carlo results. 

Consider a homogeneous Fermi mixture of two hyper- 
fine states with densities n-\ > riy. The interaction be- 
tween the two states is described by the s-wave scat- 
tering length a, while the interactions between identical 
fermions may be ignored due to the very low temperature 
of the gas. The chemical potentials \i a for the two compo- 
nents (a =t, I) are determined from n a = —dil/dfi^. We 
obtain the thermodynamic potential within the Nozieres 
and Schmitt-Rink (NSR) approximation Q by writing 
11 = ilo + Afi, where fio is the ideal gas contribution and 
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The two-body T-matrix is T 2B = 47r?i a/m, = 
e±«0 + ) is the regularized propagator for the center- 
of-mass of two atoms with total momentum K and en- 
ergy e, and (3 = l/fc^T. The single-particle proper- 
ties are described by the spectral function A a {k, u>) = 
-Im{2/[w + iO+-£fc )CT -£ CT (fe,a; + £0 + )]}, where k and w 
are the momentum and the energy of an atom in state \<j) , 
and £,k,a = k 2 /2m—[i cr is its kinetic energy measured with 
respect to its chemical potential. To include the effects of 
interactions, we calculate the selfenergy S CT (fc,u; + i0 + ) 
within the ladder approximation. The above approach 
has been shown to compare well with Monte-Carlo cal- 
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dilations for thermodynamic and single-particle quanti- 
ties 0, @, [n| , especially in the highly- imbalanced case. 

The properties of the atomic gases can be probed by 
RF photons with frequency wrf, which induce transi- 
tions from one of the hyperfme states, say |<x), to a third 
state |e) that is initially empty. Within linear response, 
the rate r <T _ >e (o;RF) is given by the imaginary part of 
the Fourier transform of the retarded spin- flip correlation 
function -i6(t-t')(&t(r,t)Mr,t),^^)A(r',t')]). 
A conserving calculation of this correlation function is a 
complicated problem [El. [Til] . However, a major simplifi- 
cation arises for the spin mixture used in recent exper- 
iments, where final-state interactions between atoms in 
states |e) and | — a) can be ignored, unless one is inter- 
ested in bound-bound transitions Then, the correla- 
tion function decomposes into Q 



dk 



A a (k, £k,a - Wrf)/(6c,ct - Wrf) 



(2) 

with /(u>) = [exp((3oj) + 1] 1 . The theory used here is 
described in detail in Refs. 0, EH . 

It is useful to introduce a measure of the fraction of \a) 
atoms bound in molecular states, which may have a finite 
or an infinite lifetime. These atoms may alternatively be 
thought of as bound in nonzero-momentum Cooper-pair 
fluctuations. For this quantity we introduce the defini- 
tion An a = —dAn/d/jia, such that n a = n° + An a with 
n° = / dk/(6;,<x)/(27r) 3 . On the BEC side of the reso- 
nance, the many-body T matrix T 2B /[l-T 2B Z(K, e+iO)} 
has a real pole at energy ek ■ The contribution from these 
real poles to An„ can easily be shown to be 



An { „ s) = 



dK Z KtCJ 
(2tt) 3 - 1 ' 
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where Zk.o — —deK /dfi a and the integral is over the val- 
ues of K for which there is a real pole. Far on the BEC 
side, An a therefore gives the density of stable molecules 
with energy ex and with an occupation probability given 
by Z K ,a- For a — > 0+, e K = K 2 /Am — l/ma 2 - (ji^ + yU T ) 
and Zk,<t = 1, such that we recover the result of n a non- 
interacting molecules. However, close to the unitarity 
limit there are no stable molecules. In this regime, the 
distinction between bound and unbound atoms is highly 
dependent on the experimental probe used to explore 
this question and the many-body state simply cannot be 
thought of as molecules mixed with atoms. 

Consider first the unpolarizcd case (n-f = n\) at 
T = 0.3T F . Here and in the following k B T F = E F = 
hp /2m = (67r 2 n^) 2 / 3 /2m is the Fermi energy of the ma- 
jority component. In Fig.[TJ we plot the fraction An a /n a 
of atoms which are bound into molecular states. In the 
BCS regime, this number is small. The fraction increases 
continuously with increasing interaction toward the BEC 
regime. At unitarity, An a ~ 0. 671,0-. On the BEC side 
of the resonance, the chemical potential becomes nega- 
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FIG. 1: Fluctuation contribution An a /n a and stable molec- 
ular fraction Any /n a for a balanced gas at T = 0.3Tp. 
Left inset: ^c/Ef vs. I/Qcfo) at T = 0.3Tf. Right inset: 
Aria/rio- as a function of temperature at unitarity. 
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FIG. 2: Spectral functions A a {k — 0, w) for a balanced gas 
at T = 0.3Tf- The frequency axis is shifted by \i a such that 
the spectral function of an ideal gas would be a delta function 
located at the origin. 



tive for \/(kpa) > 0.3, and correspondent^ the many- 
body T-matrix acquires a pole with vanishing imaginary 
part. This leads to a sharp increase in the number of 
atoms bound in stable molecular states. In fact, for 
— 0.5 we find An^ > An a , showing that, even 
if stable molecules exist, the physical picture of atoms ei- 
ther being free or bound in molecular states is too simple 
for the normal phase in the strongly- interacting limit [l3| . 
Of course, deep in the BEC regime all atoms are bound 

(s) 

in molecular states and then An a = An a = n a . 

In Fig. [2 we plot the spectral functions for the bal- 
anced gas at T = 0.37V and in Fig. [3] the corresponding 
RF spectrum for various coupling strengths. The spec- 
tral functions with k < kp, which give the main contri- 
bution to the RF signal in Eq. (J5J), look similar to the 
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FIG. 3: RF spectrum of a balanced Fermi gas at T = 0.3Tf. 
The molecular peak emerges for l/{k F a) > 0.3. The dots arc 
experimental points at unitarity and at T = 0.17> from Q. 

k = case plotted here. Deep in the BCS regime, the 
spectral function (RF spectrum) displays a single narrow 
peak shifted to negative (positive) energy, showing that 
the gas may be described as a Fermi mixture with an 
attractive mean-field shift. The spectral function in the 
unitarity limit has a double-peak structure characteris- 
tic of the pseudogap. As shown both experimentally and 
theoretically [1, [lj], there is no corresponding double- 
peak structure in the RF spectrum, which has a single 
broad peak shifted towards positive energy due to the 
effectively attractive interaction. In Fig. Owe also show 
the experimental data obtained in the unitarity limit Q . 
We see a good agreement between the theoretical and ex- 
perimental results. The discrepancy at low frequencies is 
most likely due to the theoretical curve being calculated 
at T ~ T c , while the experiment is performed at T < T c 
where the supcrfluid nature of the gas results in a true 
gap in the RF spectrum [U. From Figs. [U and f2] we see 
how the emergence of stable molecules for l/{kpa) > 0.3 
leads to the spectral function developing a gap between 
its two peaks. The peak at negative energy corresponds 
to the molecule continuum and can be understood as fol- 
lows. An atom in state \<j) with momentum k and energy 
to can pair with an atom in state | — a) with momentum p 
and kinetic energy p 2 /2m forming a molecule with energy 
Ef, + K 2 /Am. Here Eb and K = k + p are respectively 
the binding energy and the momentum of the molecule. 
Energy conservation gives 

u = E h +£ k>c - (k 7 P)2 . (4) 

In the BEC limit, A a (k,u>) therefore consists of a molec- 
ular part for w < Eb + £fc j(T separated by Eb ~ -1/ma 2 
from the atomic part. Importantly, this structure of the 
spectral function translates into a two-peak structure in 
the RF spectrum as shown explicitly in Fig. [3J The 



FIG. 4: Spectral function Aj.(fc = 0,ui) for the minority com- 
ponent of an extremely imbalanced gas at T = 0.052>. Inset: 
energy io sp and inverse lifetime l/r sp of the spin-polaron. The 
red crosses are the Monte-Carlo results from Ref. [l5( |. 

molecular part is now located at high frequencies as one 
needs to provide additional energy to break the molecule. 

In the BEC limit, the system can be considered as an 
ideal gas of atoms mixed with tightly bound molecules 
with binding energy Eb = — I /ma 2 . The RF photons 
flip the \a) atoms to the final empty state |e) with a 
rate given by the Golden Rule expression r^ e (cjRp) oc 
J dkf (^k,a)S(uJKF) / (2n) 3 . The rate coming from the \a) 
atoms bound in molecules is given by 

r^ e (WRF)0C J ^|^k| 2 n K [W(£K/2-k,- ff )] 

5(uj rf + K 2 /Am + E b - £ K /2+k, e - &c/2-k -a) (5) 

where tik is the Bosc distribution for the molecules. 
Eq. ([5]) describes molecules with momentum K splitting 
into pairs of atoms in spin states |e) and | — a) and with 
momenta K/2 ± k. The matrix element for this process 
does not depend on K due to Galilean invariance, and 
is equal to the molecule wavefunction, which is of the 
Yukawa form tp^ = V87ra 3 / (1 + k 2 a 2 ). 

The RF spectrum obtained from Eq. (O with a molec- 
ular binding energy Eb/Ep = 2, corresponding to 
l/(fci?a) = 1, is plotted with a thin line in Fig. [31 There 
is good agreement between the ideal gas model and the 
full many-body calculation for large u> since the vacuum 
Yukawa wavefunction for the molecules is accurate for 
high momenta. Closer to the threshold for molecular dis- 
sociation, the many-body RF spectrum is suppressed as 
compared to the ideal gas prediction due to Pauli block- 
ing effects; the medium suppresses the components of the 
molecule wavefunction with momenta k < kp. We con- 
clude that these interesting medium effects on the molec- 
ular wavefunction can be detected in RF spectroscopy. 

We proceed now to the analysis of a highly imbalanced 
gas (n T » ni) at T = 0.057>. In Fig. H we plot the 



4 





1 










Majority 
Minority 






1/(k F a)=-1 


















1/(k F a)=1.3 


1 

-0.3 1 0.3 0.6 0.9 
1 










Unitarity 




















-0.3 1 0.3 0.6 0.9 ' / 

1 ' , 





2 4 6 8 10 



RF offset corf^f 

FIG. 5: RF spectrum of an extremely imbalanced gas at T = 
0.05Tf. The thin lines are the ideal gas results from Eq. |(3J. 
Most of the majority spectral weight resides in the delta-like 
peak at lu-rf = 0. Only a fraction of order 1/iV (where N is the 
total number of particles) is in the high frequency shoulder. 



spectral function of the minority component for k = 
(since kp^/kp <C 1, where kpi = (Gtt 2 ^) 1 / 3 ) on the 
BCS and BEC sides of the Feshbach resonance. As for 
the balanced case, the spectral function for the minority 
component develops a broad low-energy peak for uj < Ef, 
due to the formation of molecules as described in Eq. ([!]). 
In addition, there is a low-energy narrow peak which sig- 
nals the presence of a well-defined quasiparticle, the spin 
polaron. We find that this feature is present on both 
sides of the resonance, and its energy uj sp is in very good 
agreement with recent Monte-Carlo calculations |6|, [l5| . 
Indeed, in the limit ny <C n-\ every diagram containing 
| I) holes is suppressed, and the ladder approximation 
accurately describes the main physical features of the 
system. The energy of the spin polaron is — O.GEp at 
unitarity, and on the BEC side it is pushed down by 
the molecule continuum. Spin polarons arise due to the 
dressing of minority atoms by a large number of majority 
atoms. Indeed, no trace of those is present in the major- 
ity spectral function, which consists of a narrow peak at 
lujip = due to atoms and a molecule continuum for low 
energies. In the inset of Fig. 0] we also show that the 
lifetime T sp = — l/21ml](fc = 0, lo = Lu S p) of the spin po- 
laron is very long close to unitarity. On the BCS side, 
the lifetime has a minimum (which at T = 0.05Tp lies 
at kpa ~ —1/3), and in the limit kpa — > CP it scales as 
l/(fcpa) 2 [l6j]. The lifetime decreases moving toward the 
deep BEC side since the spin polaron eventually is ab- 
sorbed in the molecule continuum and one ends up with 
an ideal gas of molecules and atoms. 

In Fig. [5] we plot the RF spectrum for the imbalanced 
case. To understand the rich physics, consider first the 
ideal gas model of atoms mixed with molecules intro- 
duced in Eq. ([5]) and given by the thin lines. The ma- 



jority spectrum is characterized by a large nonintcracting 
component, i.e., a narrow peak at zero offset, and a broad 
shoulder at wrf > 1/ma 2 , corresponding to molecule dis- 
sociation. The minority spectrum displays no signal at 
zero offset, showing that all | j) atoms are paired. Pauli 
blocking effects raise the threshold energy for molecule 
break-up in the minority spectrum by roughly Ep. In- 
deed, when a RF photon breaks a molecule and transfers 
a | J.) atom to the empty level |e), the remaining | f ) atom 
cannot occupy states in the filled Fermi sea. At T = the 
photon therefore needs to carry at least 2Ep extra energy 
to transfer the | f) atom into a state above the Fermi sea, 
but nonzero temperatures reduce this threshold shift. 

Consider now the results of the full many-body calcu- 
lation shown as thick lines in Fig. [5] In the kpa — > + 
limit, our many-body spectra fully agree with the ideal 
gas model of Eq. ([5]). Closer to the resonance on the 
BEC side, the many-body calculation agrees with the 
ideal atom- molecule mixture prediction for cjrf 3> \Eb\, 
since the large-energy uj^p 2 tails are due to two-body 
physics. However, interesting many-body effects appear 
near the molecular threshold. First, we see from the ma- 
jority spectrum that the medium increases the molecular 
energy above the two-body value —1/ma 2 . As predicted 
by the ideal gas model, the minority threshold lies above 
the majority one. Furthermore, the qualitative new effect 
of the interactions, i.e., the presence of the spin polaron, 
gives rise to a clear peak in the RF spectrum. This leads 
to the important conclusion that one can measure the en- 
ergy of this intriguing state directly by RF spectroscopy. 

We acknowledge useful discussions with A. Perali and 
C. H. Schunck, and thank B. Svistunov for kindly pro- 
viding us with their Monte-Carlo data. 
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